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We explain the reduction of the thermal conductance below the predicted universal value observed 
by Schwab et al. [Nature (London) 404, 974 (2000)] in terms of the scattering of thermal phonons 
off surface roughness using a scalar model for the elastic waves. Our analysis shows that the thermal 
conductance depends on two roughness parameters: the roughness amplitude 5 and the correlation 
length a. At sufficiently low temperatures the ratio of the conductance to the universal value 
decreases quadratically with temperature at a rate proportional to S^a. Values of 5 equal to 22% 
and a equal to about 75% of the width of the conduction pathway give a good fit to the data. 



I. INTRODUCTION 

Ballistic transport in mesoscopic system has been an active area of study. Following Landauer's approach to 
electronic conductance, several gmujps have derived expressions for the thermal conductance due to ballistic phonon 
transport in an ideal elastic beam.Eru The formula they derived shows that the only material and geometry dependence 
of the thermal conductance arises through the long wavelength cutoff frequencies of the elastic waves in the beam. 
As the temperature T — > 0, the conductance is dominated by the lowest few modes with zero cutoff frequency. The 
conductance then takes on a universal value with the value No^ir'^k'^T)/ {3h) with Nq the number of modes 
with zero cutoff frequency at long wavelengths, which is four for a free standing beam. Based on these theoretical 
predictions, great efforts have been made to observe the universal thermal conductance. Recently, Schwab et al.cl 
successfully observed the universal thermal conductance in a suspended silicon nitride bridge. Their experiment 
shows a result consistent with the universal conductance at temperatures below about 0.08 K. Above about 1 K the 
conductance rises above this value, as the modes with nonzero cutoff frequencies become excited and contribute to 
the heat transport. However in the range 0.1 0.4 K the thermal conductance unexpectedly decreases below the 
universal value. 

Motivated by the results of Schwab et al. we theoretically investigate a likely cause of the low-temperature thermal 
conductance decrease. We suggest that the conductance decrease is caused by scattering due to rough surfaces. Recent 
advanced crystal growth technology guarantees very few impurities in the material during a substrate growth, thus 
eliminating the possibility of impurity scattering. On the other hand, chemical etching can produce surface roughness 
on a scale of tens of nanometers, large enough to cause significant scattering. 

In this paper we use a simple scalar model for the elastic waves. We also use a two-dimensional approximation, 
which is accurate at low enough temperatures that modes with structure across the depth of the beam — the smallest 
dimension in the experimental geometry — are not excited. We also assume that the important roughness is on the 
sides of the beam, rather than the top and bottom surfaces, since the horizontal surfaces are MBE grown and have 
roughness at a scale of a few atomic layers, while the side faces are chemically etched. Kambili et ahlj have used a 
similar model in a numerical investigation of the effect of surface roughness on the mode propagation. 

In the following section, the details of the two-dimensional (2D) scalar model are introduced, and the scattered field 
calculation using a Green function approach is presented. In Sec. Ill, the scattering probabilities and transmission 
coefficients are calculated and the latter is incorporated into the modified Landauer formula for thermal conductance. 
In Sec. IV, the thermal conductance is evaluated numerically and compared to the experiments of Schwab et al. 

II. MODE SCATTERING 
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A. The model 



ISl 



e expression for the thermal conductance K oi a suspended mesoscopic beam connecting two thermal reservoirs 



K = 



Here the integration is over the frequency lo of the modes m propagating in the beam and lo^ is the cutoff frequency 
of the TOth mode. Also (i — l/{kBT), ks is the Boltzmann constant, T is the temperature, and h is Planck's constant. 
The effect of scattering is introduced through the transmission coefficient %n{uj): for the ideal case with no scattering 
Tjn — 1. Thus, the change of the thermal conductivity due to the rough surface is obtained by finding the transmission 
coefficient. 
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FIG. 1. 2D model used for calclation of the scattering of elastic waves by rough surfaces. 



As discussed in the introduction, we use a scalar model for the elastic waves, and model a thin geometry at low 
temperatures so that a two-dimensional calculation is adequate. Thus we consider a 2D wave-guide-like structure 
extended in the x direction and bounded at y = 0, in the absence of roughness. The waves satisfy the scalar wave 
equation, and we assume Neumann boundary conditions at the edges of the wave guide, corresponding to a stress 
free boundary condition for the elastic waves. Note that Dirichlet boundaries do not support the modes with zero 
cutoff frequency that are a crucial feature of the elastic problem. We calculate the scattering process by considering 
an elastic wave propagating in the wave guide in the +x direction with wave vector kg and entering a rough surface 
region of length L (0 < a; < L) where the rough boundaries are at y = W + f+{x) and y = f-{x) so that the roughness 
is characterized by the functions f± {x) (see Fig. Q). We assume that the top and the bottom roughness functions 
are uncorrelated, and that f± (x) is small and is differentiable. The incident wave ^'in interacts with the roughness, 
is scattered into other modes \E'sc, and leaves the rough region. The total field \E' is the sum of the incident field and 
the scattered field 

* ^ ^in + *sc. (2) 

Our task is to find an expression for '^sc and hence calculate the transmission coefficients. We do this using a Green 
function method. 



B. Green function method 

We start with the Helmholtz equation for a scalar wave at frequency u; 

VHix,y)+K^'^{x,y) = 0, (3) 
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where ^ is the total field, and K is us/c with c the wave speed. Define the Green function as a solution to the point 
sources 

V^G {x, y- x\ v') + K^G {x, y; x', y') = -6 {x ~ x') S{y - y') (4) 

with (cc', y') the source coordinates and (x, y) the observation coordinates. It is convenient to define G (x, y; x'y') such 
that it satisfies Neumann boundary conditions at the smoothed boundaries, y = O^W 

dG/dn\y^^^w = 0, (5) 

where ii is the outward-pointing normal to the surface. We then project the physical boundary conditions at the 
rough surfaces onto the smoothed boundaries to calculate the scattering. 

Multiplying Eq. (||) by G(x, y; x' , y') and Eq. (Q) by ^'(x, y) and integrating over a volume bounded by the position 
of the smoothed surfaces yields the result of Green's theorem 

^{x,y)^ f dx'^^^fi^G{x',y';x,y) 

^smooth C'n 



± / dy' 



d^{x',y') , , dG{x',y';x,y) 
G (x ,y ]x,y)--^ {x ,y) 



ox' ox' 



(6) 



where the first integral is the integration along the smoothed edges and the second is the integration over the distant 
ends taken at x — > ±cxd, and we have used the boundary condition Eq. for G to eliminate a second term in the 
first integral. 

The Green function G {x,y] x' ,y') satisfies Eq. (||). Using the completeness relation, the right hand side can be 
written 

- e^'^^^-^UkJ2My)My'), (7) 



2tt 

n 

where 0„ is the normalized transverse eigenfunction for smooth boundaries 

0n = Nn cos Xny (8) 



with Xn — rnr/W, n — 0, 1,2, and Nn the normalization factor Nn = \/2/W for n ^ and N,-,, = \/l/W for 
n = 0. The Green function is then given by Fourier transforming 

G(x yx' y') = ^ T dkY Mv)4>n{y') 

U[x,y,x,y) 27r ^ fc2 + („^/^)2 _ ^2 ' ■ W 

The k integral in Eq. is now evaluated by contour integration. The poles corresponding to propagating waves at 
k — ^ K'^ — (mr/W)^ for K > mr/W must be given infinitesimal imaginary parts ztie to yield outgoing waves. We 
then have 

G(x,,;x-,,0 = E """""?;^^^^"^^'^ (10) 

n 

where 

Using the explicit expression for the Green function the second term in Eq. (6) can be shown to be just the incoming 
wave (x,y), so that 

^se{x,y)^ f dx'^^)fi^G{x',y';x,y) (12) 

J smooth 
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C. Boundary perturbation 



In the absence of roughness the field 4" satisfies Neumann boundary conditions at the smooth boundary, and so the 
scattered field would be identically zero as expected. Correspondingly, for a rough surface with small f±{x) we can. 
calculate d'^/dn at the smoothed surface appearing in the integral by expanding about the stress- free rough surface.H 
We will present the calculation for the rough lower surface, and simply double the scattering probabilities assuming 
uncorrelated roughness on the two surfaces. 

Firstly, express the unit normal vector as 

n= -y + f'_ (x) X. (13) 
Then impose the Neumann boundary condition at y = /_ {x) 



d-^{x,v) , {x,v) 



= 0. (14) 

y=f- (^) 



Now expand this equation about y = in terms of /_ and retain only terms that are first order in / and /'. This 
gives for the normal derivative at the smooth surface up to first order in /, /' 

5„* (a^,y)l,=o = i^^v) - f-dl^ i=^^y)] \y=o ■ (15) 

Thus the scattered field to first order in the roughness amplitude is 

(x, y) ~ y dx'G {x', y'; x, y) [-/„ {x') a^,*;^ {x' , y') + f__ {x') a,,*.^ (a;', y')] |^,^^ (16) 

where we can replace the field appearing in the integral by the incident field at this order. 

It is now straightforward to insert the explicit expression for the Green function Eq. ( [l0|) to calculate the scattering 
from a normafized incident wave entering in the m th mode 'I'in {x, y) = 5'm(2;, y) = Nm cos [xmV) e*'^"'^: 

vf^c {x, y)^ f dx' ^ ^^^^ cos ixny) eM^--'\ [/_ (x') xl + ^k^f- {x')] e^^-^' . (17) 

D. Scattered field 

Outside the scattering region, we may take an asymptotic form for the scattered field Eq. ([l7| ) 

M'sc(a;-+oo,y) = 5^e''="-cos(x„y) / dx' [xlf- {x') + ik„jL {x')] e^'^^--^-^^' , 



d^'^^t^ [xlf- i^') + ^kmfL {x')] e^('='"+'='.)- (18) 

n °° 

giving the forward scattered field and back scattered fields respectively. The terms in can be simplified by 
integration by parts 

dx if^ {x ) kjYiE^ mT 'n) — (/cin ^ kji) /— {k„i ^ fc„) , (19) 

where /_ is the Fourier transform of and we have used the fact that the roughness is confined to < x < L so 

that /-(±oo) = 0. 

Now using K = uj/c = ^ Xm + J where c is the velocity of the elastic wave, we get 

^sc {x ^ ±oo, y) = ^ ^^^^/- ((fc™ T kj) {K^ T k^k^) {x, y) , (20) 

n 

for the forward and backward scattered waves, expressed as a sum over normalized waves 
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III. THERMAL CONDUCTANCE 



Let t±n,rn be the scattering amplitude from mode m to ±rt, where the plus sign is for forward scattering and the 
minus sign is for back scattering. Then 

t±n,m — ^TT (^K ^ fen^m) /— {km "F ^n) • (21) 

To calculate the transmission coefficient appearing in the expression for the thermal conductance we need the energy 
flux scattering probabilities a±n,m given by multiplying |t„,mp by the ratio of the group velocities 

<7±n,m = T~ (\t±n,mf^ , (22) 

where we can also now average over the ensemble of surface roughness represented by the angular brackets. This 
finally gives 

a±„„ = [K^ T knkm] S^~gikm T kn)L (23) 



using 



f^{k)^\ ^ 5'~g{k)L, (24) 



where 5'^g{k) is the Fourier transform of the surface roughness correlation function with S the roughness amplitude. 
For the characterization of the rough surface, we assume a Gaussian correlation function g{x) = 1°" where a is 
the correlation length of the roughness, so that 

~g{h) = V^aexp [-a^P/A] . (25) 

To calculate the thermal conductance we must recognize that not all scattering processes decrease the heat transport. 
A wave entering in mode m has four possible outcomes: after the scattering events it may stay in mode m propagating 
forward; it may be converted to mode n also propagating forward; it may stay in mode m but propagating backward; 
and finally it may be converted to mode n and propagating backward. The former two cases do not change the 
heat transport, since each mode at frequency u contributes the same amount to the conductance. The two back 
scattering events do reduce the heat transport however. Thus the backward scattering rate a-n,m contributes to the 
reduction of the thermal conductance while a+n,m is the coefficient for forward scattering, and leaves the conductance 
unchanged. We therefore define the conductance attenuation coefficient per unit length of the rough surface waveguide 
as 7m = (2/i) J2n '^-n,m (whcrc the factor of 2 is to include the scattering off the top surface) 



n 



The conductance attenuation coefficient 7™ gives the exponential decay rate of the wave in mode m, so that over 
a length L the transmission is 

T„ = e--^-^. (27) 
To calculate the thermal conductance at a given temperature, we insert Eqs. (pq),d2^) into Eq. (nl): 



1 .00 ^2g/3R^ 
^ = y .T^Z dajJTTm:, TT7 c^p 



kgT^ ^ 27r 7^,„ {el^^'^ - 1)^ 



(28) 



The contributions to the conductance attenuation coefficients per unit length 7,„ for the first few modes are shown 
as a function of the mode frequency in Fig. |^. A transverse correlation length of a/W — 0.75 was used in the figure. 
The backscattering amplitude from the lowest mode (mode 0) to its reverse is 

7oo(-,a,^) = 2.V2^^,-.^.Vc^ (29) 
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This expression is finite for all frequencies. It has a maximum at a frequency lo = c/a depending on the roughness 
correlation length, with a peak value of order [5"^ /aW"^). The higher modes have a diverging back scattering propor- 
tional to {lo — ujm)^^ at the cutoff frequencies ujrn- In addition each 7^ has a contribution diverging as {lu — a;„)~^/2 
at the onset of the rtth mode. These divergence are due to the flat spectrum at the mode cutoff frequencies, and will 
also be found in a full elastic wave calculation. 
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Reduced Frequency co/A 

FIG. 2. Scaled attenuation coefficient {VV^ /S^aL)\a-n,m\'^ as a function of reduced frequency, lj/A where A = -kc/W with 
c the average velocity of the elastic waves: solid - from mode m = to mode — n, n = 0..3 and dashed - mode m to mode — m, 
m = 1..3. A value of the roughness correlation length a = 0.75W was used. 



At low enough temperatures only the lowest mode with fco = lo/c contributes to the thermal conductance, and 
only the backscattering of this mode given by Eq. ( ^9| ) reduces the conductance below the universal value. This 
reduction is plotted as a function of the temperature scaled by hc/kBO, in Fig. ^. The temperature of maximum 
reduction depends on the roughness correlation length a, whereas both the roughness amplitude and correlation 
length change the magnitude of the reduction. For small roughness, we can expand the exponential term in Eq. ( [2^ ) 
and at low temperatures only small w contributes to the integral so that the Gaussian factor may be replaced by 
unity, exp[ — a^w^/c^] ~ 1. This leads to 



K 



3h 



8^9/2 S^aL fkuT 



1^4 V /lA 



(30) 



where A = -kc/W is the spacing between the mode cutoff frequencies. Thus at low temperatures, the conductance di- 
vided by temperature should show a quadratic temperature decrease with an amplitude depending on the combination 
of roughness parameters 5^a/W^. 



IV. COMPARISON WITH EXPERIMENT 

To compare with the experiments of Schwab et al.i we use the following geometry and material parameters. We 
take a wave guide structure of rectangular cross section with width W — 160 nm and length L = \ fjjn.. In the 
experimental geometry the width varied along the length to provide smooth junctions with the reservoirs. This was 
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done to eliminate scattering off abrupt changes in the geometry. We use the width at the narrowest point as our 
estimate. For the length we use the length of the central portion over which the width is fairly constant. Since the 
length only occurs in the combination S^L, changing the value of L used will only change the estimated value of 6. 
We use a wave propagation speed c — 8250 m/s which is the average of the velocity of longitudinal and transverse 
elastic waves in silicon nitride. 

The roughness parameters are not known a priori. As a first attempt we might try to estimate the combination 
S'^aL/W^ from the quadratic decrease in the thermal conductance at low temperatures, Eq. (|30|). This would give the 
value aS'^L/W'^ ~ 0.05. However, from Fig. ||we see that the quadratic low-temperature fit is only good up to about 
a quarter of the temperature of the maximum backscattering of the first mode. If we estimate this temperature from 
the minimum in the measured conductance, we find that the data does not extend to low enough temperatures to 
provide a reliable fit, and so this value can only be used as an order of magnitude. In fact our "best fit" (see below) 
over temperatures up to 1 K corresponds to a value aS^L/W^ about a factor of 4 larger 

It is interesting to use this value of the roughness parameters to estimate the strength of the scattering of the 
higher modes. For example, for the first mode, with cutoff frequency A, and at a wave vector n/W corresponding to 
a frequency ^/2A we find for the backscattering into the same mode 

7i(fci = n/W)L - 16exp {-tt^o^/W^) . (31) 

The scattering increases for smaller wave vectors, diverging at onset as shown in Fig. ^ Remember that the trans- 
mission amplitude is e""*"^^. This means that the scattering of the higher modes is strong over the 1^ length, unless 
sufficiently reduced by the exponential factor arising from the reduced roughness at short length scales. To fit the 
higher temperature data using Eq. (|28| ) we will find that we need a value of a comparable to W. Although this 
strongly reduces the value of 7i(fci = t:/W), there remain frequency ranges where the scattering of this mode and 
other modes is strong. An interesting consequence is that a significant fraction of the thermally excited phonons at 
temperatures of order 1 K are predicted to be localized in the experiments of Schwab et al., with a localization length 
less than the length of the bridge. Unfortunately, in this regime the estimate of the contribution to the conductance 
from these modes predicted by our lowest order scattering calculation, will not be accurate. 




Scaled Temperature T 

FIG. 3. Reduction in the thermal conductance divided by temperature due to back scattering of the lowest mode, expressed 
as the ratio to the universal conductance divided by temperature and then scaled by aW^/S^L, as a function of temperature 
scaled by hc/kBa. 

From Fig. ^ we can suggest two mechanisms that might account for the observed minimum in the dependence of 
K/T on temperature. The first mechanism ascribes the minimum in K/T to the behavior of the first mode alone, as 
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plotted in Fig. The upturn in K/T arises from the reduced scattering of the lowest mode as the wave vectors of 
the important modes increase with temperature. The second mechanism supposes that the scattering of the lowest 
mode is responsible for the decreasing K/T at low temperatures, but that the subsequent increase is from the thermal 
excitation of the higher modes. For our "best fit" values of a, 5 (see below) the results are summarized in Fig. IJ. 
The picture is quite complicated, with both the reduced scattering of the lowest mode and the thermal excitation 
of the higher modes contributing to the rise in K/T with increasing temperature. Furthermore, due to the strong 
scattering of the higher modes near their cutoff frequencies, these modes become important in the transport at a 
higher temperature than would be estimated simply from their cutoff frequencies. The higher modes excited near 
their threshold frequencies are localized and do not contribute significantly to the transport. 




Scaled Temperature 

FIG. 4. Contributions to the thermal conductance K divided by the universal value Ku from the first few modes for the 
ideal no scattering case, and for the rough case with scattering as a function of the scaled temperature fcsT/fiA: solid line 
- total thermal conductance for the rough surface case; dotted line and short-dotted line - conductance of mode (ideal and 
rough); dashed line and short-dashed line - conductance of mode 1 (ideal and rough); dash-dotted line and short-dash-dotted 
line - conductance of mode 2 (ideal and rough). Values of the roughness parameters were a/W = 0.75 and 5/W — 0.22. 

In Fig. ^ the thermal conductance calculated using Eq. ( |2^ ) is plotted together with the ideal (no-scattering) 
conductance and the measurements of Schwab et al. The conductance is scaled such that the universal conductance 
appears as unity. The roughness parameters a/W = 0.75 and 5/W — 0.22 (so that a5'^L/W^ — 0.23) were used, 
and yield a reasonable fit to the data. Our 2D model shows the same trend as the experimental data: a decrease 
in the thermal conductance below the universal value at low temperatures where only the lowest modes are excited, 
then a gradually increasing conductance as other modes are excited and the scattering of the lowest mode is reduced. 
Comparison to the ideal (non-scattering) curve shows that the scattering is important over the whole temperature 
range examined (T < IK). These values of (5 = 35 nm and a = 120 nm appear reasonable when one considers 
the physical process of constructing the mesoscopic bridge structure. For example, a typical chemical etch of silicon 
nitride can easily produce a few tens of nm in roughness amplitude. Electron micrographs of the actual structure 
used in the expcrimentQ show roughness on scales comparable to the ones we estimate. 

There are small but systematic differences at very low temperatures, where the conductance is dominated by the 
lowest modes, and the theory should be most accurate. The discrepancy suggests that we are overestimating the 
scattering at long wavelengths. A roughness spectrum g(k) ^ fc^e^" * with a reduced amplitude at small wave 
numbers gives a better fit to the data. Such a form might be physically reasonable, since we might expect the 
roughness to be largest at a scale of order the minimum dimension of the structure, and reduced at larger scales than 
this. However, since the scalar model does not account for the mode structure of the elastic beam accurately, it is 
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probably unwise to use the discrepancies in Fig. ^ to make any firm deductions. Such conclusions must await a more 
accurate treatment of the modes within elasticity theory. 




Temperature T (K) 

FIG. 5. Thermal conductance relative to the universal value Ku as a function of temperature for the ideal case (dotted line), 
the rough surface case (solid line), and the data of Schwab et al. (circles). The roughness parameters used were a/W — 0.75, 
S/W = 0.22. 



V. CONCLUSION 

We have investigated the cause of the thermal conductance decrease below the universal value at low temperatures 
by employing a Green function approach to calculate the reduced transmission of the elastic waves due to surface 
roughness, and then using Landauer's formula for the thermal conductance. 

At low temperatures, the conductance divided by the temperature is dominated by the lowest mode. The scattering 
of this mode reduces the conductance below the universal value with a quadratic dependence on temperature for low 
temperatures with an amplitude proportional to the combination of roughness parameters aS^. As the temperature 
increases, higher modes begin to play a role and the scattering of the lowest modes is reduced, so that the conductance 
increases. We find that the effect of scattering is always significant, reducing the conductance below the ideal ballistic 
value over the whole temperature range we investigate T < 1 K. Considering the simplicity of our model our results 
agree well with the experiment of Schwab et al. In future work we will present results for a full elastic theory treatment 
of the thin bridge. 

ACKNOWLEDGMENTS 

The authors are grateful to Keith Schwab for providing the data and an electron micrograph of the experimental 
structure, and Miles Blcncowc for carefully reading the manuscript and providing useful suggestions. This work was 
supported by NSF grant no. DMR-9873573. 



9 



^ D.E. Angolcscu M.C, Cross, and M.L. Roukes Superlattice and Microstructure, 23, 673 (1998). 
^ L.C. Rego and G. Kirczonow, Phys. Rev. Lett. 81, 232 (1998). 
3 M.R Blencowc, Phys. Rev. B 59, 4992 (1999). 

" K. Schwab, E. A. Henriksen, J. M. Worlock, M. L. Roukes, Nature 404, 974 (2000). 
5 A. Kambih, G. Fagas, V. I. Fal'ko, and G. J. Lambert, Phys. Rev. B. 60, 15593 (1999). 

^ S.V. Biryukov, Yu.B. Gulyaev, V.V. Krylov, V.P. Plessky, Surface Acoustic Waves in Inhomogeneous Media, Springer- Verlag, 

BerUn (1995). 
^ K. Schwab, Private communication 



10 



